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Summary 

A solution is presented for the problem of a finite-length crack branching off the interface between 
two bonded dissimilar isotropic materials. Results are presented in terms of the ratio of the energy release 
rate of a branched interface crack to the energy release rate of a straight interface crack with the same total 
length. It is found that this ratio reaches a maximum when the interface crack branches into the softer 
material. Longer branches tend to have smaller maximum energy release rate ratio angles indicating that all 
else being equal, a branch crack will tend to turn back parallel to the interface as it grows. 


1. Introduction 

The interface crack problem was addressed in the late 1950 s by Williams ( 1959). Williams used an 
eigenfunction expansion technique to solve the interface crack problem and discovered complex singularities 
and rapidly oscillating stresses near the crack tip. Williams also estimated the oscillatory region to be very 
small. 

In the mid 1960’s, various workers [England ( 1965), Erdogan (1965), Rice and Sih (1965)] came up 
with closed form solutions to the interface problem. These solutions verified Williams' (1959) finding of 
complex singularities and rapidly oscillating stresses near the crack tip. England (1965) also pointed out 
that the complex singularities lead to the physical impossibility of the crack faces interpenetrating near the 
crack tips. Like Williams (1959), these workers found the region where these phenomena occur was very 
small. 

In the mid 1970’s Comninou ( 1977) put forth an alternate solution that addressed the difficulties at tire 
crack tip. Comninou resolved the problems at the crack tip by solving tire interface crack problem 
assuming there is a small contact zone near the crack tips. This zone was found to be extremely small. 
While calculations were carried out for a limited range of Dundurs’ constants (Dundurs 1969), Comninous 
work provides a rational explanation of the oddities occurring near the crack tips. 

Lately, the interface crack problem has received renewed attention. Gautesen and Dundurs (1987) have 
solved the problem formulated by Comninou (1977) using quickly convergent series. Unlike 
Comninou’s (Comninou 1977) results which were carried out for only severe material mismatches, 
Gautesen and Dundurs obtain results covering all material combinations. Symmington (1987) has 
completed Williams' (Williams 1959) eigenfunction analysis by finding a set of integer eigenvalues left out 
by Williams. Rice (1988) has proposed that the peculiarities at the crack tip are not too critical. In an 
argument similar to the concept of small scale yielding, he asserts that the complex stress intensity factors 
will be indicative of the general state of the crack tip even if they do not correctly represent the state 
immediately surrounding the crack tip. Rice also provides a framework to interpret the complex intensity 
factor in terms of the classical form stress intensity factors K\ and Kn- He also includes an interesting 
discussion of the interaction of the crack length and load phase angle interaction. Hutchinson, Mear, and 
Rice ( 1 987) have looked at a crack completely in one body but very near the interface. Their solution, 
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based on dimensional arguments, energy considerations, and the representation of tlie crack as a distribution 
of dislocations, examines criteria for crack growth parallel to the material interface. Park and 
Eamtme ( 1 986) use various integrals to describe the interface crack tip characteristics and include a 
discussion of the properties of these integrals. Shih and Asaro (1988) have investigated the interface crack 
with plasticity. They find that the pathological features of the linear elastic analyses are lessened by non- 
linear behavior and that the linear elastic solutions are good in regioas not immediately around the crack tip. 
Delate and Erdogan (1988) have looked at the interface crack by modelling the interface as a transition 
material; they put the crack in a thin third layer between the two half-spaces. They present some results 
regarding the direction of maximum Kj of the crack. Generally, Kj is larger as the crack tends toward the 
softer material. He and Hutchinson (1988) have performed an analysis very similar to the present one for 
the case of a branched semi-infinite crack. They have presented extensive results on the initial branching 
angle, couched in terms of strain energy release rates. 

In this paper, the branching of a finite length crack is analyzed. For small branch lengths the 
asymptotic results of He and Hutchinson are reproduced. Consideration is also given to further branch 
growth and the associated influence of a more global portion of the stress field. Based on the obtained 
results, some conjectures are made concerning the fate of branching cracks. 

2. Formulation 

The problem configuration is depicted in Figure 1, which shows the main crack, the branch crack, the 
loading at infinity, and the coordinate system. In the upper half-plane (S | ), the shear modulus is p i and 
Poisson s ratio is Vj. In the lower half-plane (S 2 ), the shear modulus is P 2 and Poisson’s ratio is v 2 . The 
boundary conditions for this problem require that both the main and branch cracks are traction free and that 
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at infinity the stresses approach a yy and T yx . Single-vaiuedness of the displacements at infinity is also 
required 

Some techniques useful for dealing with interface boundary conditions which are used in the problem 
formulation are presented first. These results follow directly from the formulation used by Clements (1971) 
to treat interface cracks in anisotropic materials. 

In terms of the Muskhelishvili (1953) potentials, the stresses and displacements in isotropic elastic 
bodies may be expressed as follows: 


(Oyy - ityx)i = <I>i(z) + <Di(z) + Z <J>;’(z) + 'Fi(z) (1) 

(tfyy + ^xxli = ^[^(z) + (z)j (2) 

2 Mil-- - i - = KiOi(z) - + zO^) + 4^(3 (3) 

\dx / * 

in which the subscript i (i = 1,2) denotes "in region i"; correspond to the potentials for the upper 

half plane and <J> 2 , ^2 correspond to the potentials for the lower half plane. Moreover, z Is the complex 
variable x + i y, the prime denotes differentiation with respect to z, an oveibar denotes conjugation, 

3 - v 

k : = 3 - 4v for plane strain, and k = -for plane stress. 

1 + v 
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For interface problems, it is convenient to introduce add itional potentials as follows. Making use of 
the fact that if f (z) is analytic for z in region R, then f (z) = f (z) is analytic for z in region R, the 
following analytic jump potentials are constructed: 


^ = l & 1 (z) - [<D 2 (z) + z Q 2 ' (z) + »P 2 (z)1 z g Si 

1 <& 2 (z) - ki (z) + z<Di (z) + 'Fi (z)j z e S 2 


(4a) 


Oq = 


I ^-®1 (z) +^-[<I> 2 (z) + z4>2 (z) + *P 2 (z)j 

_/ 2 Pi 2(i 2 

|-^<I» 2 (z) + zO,*(z) + l (z)] 

\1 p 2 2 (0,1 

For interface problems, the interface conditions are usually of the form: 


Jump in Displacement = 


z g Si 
z g S 2 


(4b) 


~ 1 ^yxfe 

- (Oyy 

* T yx)l 



(5a) 


f 3u 

. dv\ 

/du 

,dv| 

*] 


(5b) 


,dx 

u 
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Note that in terms of these jump potentials, the interface conditions can be expressed as: 

Jump in Stress = Os2 (x) - Osi (x) 

Jump in Displacement = Qd 2 (x) - Odi(x) 

The inverse relationship between the jump and Muskhelishvili potentials is given by: 


(6a) 

(6b) 


Ol(Z) = Qj 
'Fl(Z) = Q2 
3MZ) = Q 2 
4^2 (Z) = Qj 


2p 2 


Osi(Z) + Odi(Z) 


Os 2 (Z) + Quiz) 

2 |I2 


Q S2 (Z) + Qd2(Z) 


<D,(Z) - z4>,(Z) 


2 pi 


V K ' Osi(Z| + Q D,(Z| 

2 Pi 


^(Z) - z02(Z) 


(7a) 

(7b) 

(7c) 

(7d) 
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where 


Q = 2 ^2 

° + Ki\L 2 

(h = _2^__ 
^2 + *2^1 


(8a) 


(8b) 


This formulation allows for straightforward development of the various interaction potentials derived below. 


3. Solution 

The solution is obtained by a Green’s function technique based on distributing dislocation singularities 
along die branch. The solution for the interaction between an interface crack and a dislocation is found by 
superposing the solutions for: (1) a dislocation in S\ near an interface; (2) two perfectly bonded semi- 

infinite bodies loaded at infinity with Oyy and T yx ; and (3) an interface crack loaded with the negative of the 
stresses produced by (1) and (2). This solution for a single dislocation interacting with an interface crack 
(with the appropriate far-field loading) is then modified to model a distribution of dislocations along the 
branch, and by requiring zero stress along the branch face a singular integral equation is obtained for the 
unknown dislocation distribution. The relevant physical quantities can be calculated after numerically 
solving the integral equation. The next three sections present the solutions to each of the problems 
described above. 


3.1. Dislocation Near an Interface 

The solution to this problem is well-known (e.g. Head, 1953, Dundurs and Sendeckyj, 1965) , but it is 
convenient to re-derive it here using the jump potential approach. Consider the problem of a dislocation 
located at point zq in the upper half plane, which is assumed to be perfectly bonded to the lower half plane. 
The interface boundary conditions are: 


or 


(Gyy 

■ iT yx)i 

- (Gyy 

" i ^yx 


.dv\ 



.dv' 
— 1 

tdx 

‘aJi 


[dx 

dxj 

Osi 

■ oS = 

o 1 

! l x l 

< », 

Qdi 

Ore = 

0] 

1 I I 



2 |x| < oo t y = 0 


(9a) 


(9b) 


The solution is sought in the form: 

- D » D - D 

^1 “ ^1, Singular + ***1, Continuation 

^1 ~ ^hSingular + ^1, Continuation 


(10a) 

(10b) 
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Ct>2 — *1*2, Continuation (K)c) 

u/ D - u/ D 

*2 — *2, Continuation (lOd) 

in which ^singular and ^Singular are the (known) potentials for a dislocation in a homogeneous full plane, 

and ^Continuation and ^Continuation are additional (unknown) potentials needed to satisfy the continuity 
conditions along the interface. Starting off with the full plane dislocation solution in region 1 and all 


region 2 potentials equal zero: 


_ A 

^1, Singular ~ z - Zq 

( 1 la) 

_ A i a ^0 

r 1, Singular - z - Zo + A ^ ^ 

(lib) 

^*2, Singular ” ^2, Singular = 0 

(11c) 


where A = p { e i0 { [u r ] + i [vq] }/irc(Ki + 1), where [u r ] and [v©] are the jumps in the tangential and 
normal displacements accrued upon circling the dislocation. 

The corresponding il potentials are: 


_ A 

**S1, Singular — z - Zq 




^S2, Singular — 


Z - 20 


^=r + A 


ZO - ZQ 

(z - z of_ 


(12a) 

(12b) 


qD, 0 . , _ *i A 

^A)!, Singular — 2 Z - Zq 


(12c) 


_ 1 


Op2. Singular — ^ 


Hi 


Z - 70 


+ A 


zq - zp 

(z - 


(12d) 


Obviously, the boundary conditions given by equation (9b) are not satisfied. To patch up the interface 
boundary conditions, the following potentials are added, which are analytic in their respective regions: 


^^1, Continuation 

~ ^S2, Singular 

(13a) 

^S2, Continuation 

= ^Sl, Singular 

(13b) 

^Dl, Continuation 

= £^D2, Singular 

(13c) 

^D2, Continuation 

= ^Dl, Singular 

(13d) 


Using equation (7) to invert these to Muskhelishvili (1953) potentials gives: 
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(14a) 


<D 


D 

1, Continuation 


= A 


(a - p) 


.d 


^1, Continuation — A 


Continuation — A 


^2, Continuation — A 


[(i + p)(z - 

■ zo). 

(a + p) 

1 

(l - p)(z - 

zo). 

1 + a 


(l - p)(z - 

zo) 

[(1 + a)(z5 

- Zq) 


4* A 

(a - p)(zo - z«) 

T /A 

(l + p)(z - zo) 2 


. D - D» 

^1, Continuation “ z Continuation 


+ A 


1 + a 


1 + p)(z - zoPJ l(l + p}(z - zo) 


- D D* 

t 2, Continuation “ z ^2, Continuation 


(14b) 

(14c) 


(14d) 


where a and p are Dundurs 1 (Dundurs 1969) constants defined by: 


_ M-2 (1^1 + 1) - M-i (tc 2 + l) 
a m- 2 (ki + 1) + M-l(K2 + 1) 

(15a) 

_ |1 2 (K) - 1) - |Xi (k 2 - 1) 

M- 2 ( k i + 1) + M-i ( k 2 + 1) 

(15b) 


3.2. Bonded Semi-infinite Half Planes Under Remote Stress 

The boundary conditions for the interface body (no crack) under remote stress are: 

12s i - £2<$2 — 0 \ | | n 

M ! I x I < °°, y = 0 ( 16 ) 

- &D2 = 0/ 

As pointed out earlier, £2 S 1 , £2 D | and £2 S 2, &D2 are analytic in their respective regions. Moreover, since 
these two potentials are equal along the region boundary (as given by the above boundary condition), it can 
be concluded that £2 S and £2 D are analytic everywhere. By Liouville’s theorem, since £2$ and £2^ are analytic 
everywhere and bounded, they must be constant. Thus, set: 


Os = Cs = RCs + i ICs 

(17a) 

£2d = Cd = RCd + i ICd 

(17b) 


in which RCs, ICs, RCd, and ICd are real constants. 

Inverting relations (17) to Muskhelishvili (1953) potentials gives: 
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0>i = Q! 


2^ CS + CD 


Vi = Qi 


-k 2 

2\i 2 


Cs + Cd 


- Qi 


4 Cs ' Cd 


®2 = 02 


2 (Xi 


Cs + Cd 


¥2 = Qi 


-Ki 


2 |Xi 


Cs + Cd 


- Q 2 


2^ Cs + Cd 


The two complex constants Cs and Cd are determined from the conditions at infinity 


o~ - it- = or + o>r + *r 

^ + <Ci = 2 Ur + or) 

+ <^2 = 2 ( ) 


Inverting these leads to the following expressions for the potentials 


oo 00 

a yv + 

O] = - ^ — - — + 1 

4 

OO oo 

” " nl + it; 

1 n 


<P - j)$ + J_^x 


ICs 


y* 


oo Gyv + Go 
02 = " - — + 1 


(ft + ilv + 1 + a 


ICs 


Gyy * G^O 

T 2 = — - — + IT, 

2 


y* 


(1 8a) 
(18b) 


(18c) 


(18d) 


(19a) 

(19b) 

(19c) 


(20a) 

(20b) 

(20c) 

(20(1) 


in which the arbitrary constant 1C$ does not affect the stresses and is related to a rigid body rotation. For 

a - B 

convenience, it is chosen so the lmag(<I>| ; Far pj e |d) = 0. This gives lmag(d> 2 , p ar Field) =~ x yx . 

In determining the relations (20) it is also found that the x-direction stresses are related: 


4 P - 2 a ^ 1 + a 

** = ~T~r c T a >y + TT^ 


(21) 


This result for g xx2 and o xx \ is typical of interface problems (see e.g. Dundurs (1969)). Since this analysis 
will be coasidering branches which are not strictly parallel to the x-axis, these x-direction stresses are 
important: they can not be ignored as in the normal interface crack problem. Since it is generally not 
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possible to have these stresses vanish in both regions simultaneously, some choice must be made for their 
values. The choice made here is to take o xx 2 = ~ cr xx j , which implies that the net force on a vertical cut is 
zero. It was verified using a finite element analysis of bonded interface bodies of finite dimensions that 
there exists a region at the center of the specimens near the interface where the tensile stress parallel to the 
interface in one body is equal and opposite to the tensile stress parallel to the interface in the other body. 
Thus the present choice of x-direction stresses would be appropriate for interface cracks far from free 
boundaries with no applied x-direction tractions. 

3.3. Interaction with the Main Crack 

The stresses due to the dislocation near an interface and the far field loading are removed from the crack 
face by the following potentials: 

C C C 

^ = ^Dislocation + Field (22a) 

C C C 

^ = ^Dislocation + ^Far Field (22b) 

C C 

where location amoves the stresses on the crack face due to the dislocation and d> Far FjcJd removes the 

stresses due to the far field loading. To obtain these potentials, first consider the interface crack under 
arbitrary crack face loading. The standard solution has been presented by many investigators, e.g. England 
(1965), Erdogan (1965), and Rice and Sih (1965), but for convenience will be rederived here via the Jump 
potentials in a manner similar to Clements (1971). 

For a crack on an interface, the interface boundary conditions are: 


(Oyy - 

O 

ll 

£ 

1 

£ 

1 

| x | < 00 , y = 0 

(23a) 

(u - 

i v')i - (u - i v'fe = 0 

|x| > C 

(23b) 

Oyy - 

i "tyx = AM 

|x| < C 

(23c) 


These boundary conditions in terms of the standard Muskhelishvili (1953) potentials are somewhat 
cumbersome. However, as shown by equation (6), in terms of the jump potentials the first boundary 
condition is simply: 


Osi(x) - Q s2 (x) = 0 |x| < «. (24) 

Using an argument similar to the one in the Far-Field solution, since £ 1 $ is analytic everywhere and 
bounded, by Liouville's theorem it must be a constant. Moreover, for zero stress at infinity, f2 s = 0. 

The second boundary condition in terms of the jump potentials is: 

£2 D i(x) - Od 2 (x) = 0 |x|>c (25) 

The thir d boundary condition in terms of the jump potentials is: 

QlO D1 (x+) + CkO^x-) = f(x+) = f(x) (26a) 
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02^1>2(x ) + Ql^Dl(x + ) = f(x ) = f(x) 


(26b) 


or 


£2d + mO D = f(x) 

Qi 

_ _ Q2 _ 1 + P 

m = — = 

Qi 1 - p 

These last two boundary conditions (equations 
known solution: 


(27a) 

(27b) 

and 27a) define a Hilbert problem, with the well- 


Q d ( z ) 


_L 

Q 


f(x) 


Qi X + (x)(x - z) 


dx + X(z) P(z) 


(28) 


in which 


X(z) = (z - c)^ ' 1 (z + c ¥* 


Y= i._L logm 

= 1 + i £ 


C = 2 k ‘° g - 


and P(z) is a polynomial determined by far-field behavior. 

After is obtained, the jump potentials can be inverted back to standard potentials using equation (7). 
In the following two subsections the general result in eq (28) is specialized for the dislocation interaction 
and the far-field loading. 


3.3.1 Removal of Dislocation Stresses from Crack 
To remove the stresses on the crack caused by the dislocation solution, equal and opposite tractions are 
added. The stresses due to a dislocation near the interface are obtained by substituting z = z = x into (1). 
These stresses are: 


<*yy ” i ^yx — A 


+ A 


1 + a( 1 

1 dU-Zq 


1 + a 


1 - P ^ 1 + P' x 

1 + a (zo ~ zd) 

1 + p (x - 


y 


(29) 
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Substituting (29) with a negative sign into (28) and performing the necessary integrals enables the 
interaction £2 potentials to be determined: 


^S, Dislocation — 0 


(30a) 


0 C _ - 2 A 1 

^“D, Dislocation ~ { + m 

- . 2 ,^ — ^-(zo - 
1 + m Qi 

+ X (z) P(z) 


i-+-&F(z,zo) + J-i«F(z,zo) 

1 - P 1 + P 

1 + — G (z , z5) 


1 + P 


(30b) 


where 


F (z , a) = - 


1 


G (z , a) = 


2 (z - a) 

dF (z , a) 
da 


1 - 


X (z) 


X (a) 


(31a) 


(31b) 


Inverting these to standard potentials gives: 


, Dislocation 


, Dislocation “* 


2 A [ 

1 + fi F(z,zo) + 1 +SL 

1 + m 

1 - P 

1 + P 

2 A 

(zo - ^)) 

* + a G (z, z5) 

1 + m 


.1 + P J 

2mA 

(zo - zo) 

- -- g G (z, 

1 + m 


1 + P 

2mA 

1 + a F (z, zoi + 1 ± - g 

1 + m 

1 ' P 

1 + P 


F (z, zq) 


F (z, zo) 


_ c . c • 

9^1, Dislocation " z ™1, Dislocation 


(32a) 


(32b) 


<t>2 , Dislocation 


2mA 
1 + m 


1 + g F (z, zo) + 
1 - P 


2 mA ^ . zjj) 
1 + m 


1 + ~ — F (z, zo) 

1 + P 

1 + - G (z, z5) 

1 + P 


(32c) 
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^2, Dislocation ~ ~ ^ ^ (A) ' Zo) 

1 + m 


2 A 


1 + g 

1 - P 

c c ' 

- Didnrahnn ” zO?n 


1 + in 


1 ^-G(z,zo) 

1 + P 

F (z, zo) + 1 + (x F (z, zo) 

1 + P 


(32d) 




Each of the potentials in (32) contains an X(z) P(z) term following inveision . However, requiring that 
the stresses vanish at infinity and requiring that there be no net displacement at infinity (no net 1/z term) 
leads to all the P(z)'s being zero. 


3.3.2. Removal of Far Field Stresses from Crack 

This portion of the problem is identical to the classical interface crack problem. Proceeding as usual the 

oo oo 

crack is loaded with - (a yy - i x yx ) such that the net stress on the crack is zero when this solution is 

superposed with the remote stress solutioa After performing the necessary integrals in equation (28), the 
potentials are: 


^S, Far Field = 0 


^D, Far Held 


yy 


it, 


y* 


' yy 


* "^yx 


(1 + m)Q! 
P(z)X(z) 


(1 + m)Q! 


[z + (2y - l)c]X(z) 


(33a) 

(33b) 


Again, P is deteimined by far field behavior and turns out to be zero. Inverting to standard potentials 
gives the following result which is of course identical to the previous solutions described earlier: 


. C 

^1, Far field 


'F 1, Far Held 


< t*2. Far Reid = 


C 

'F 2, Far Reid 


V-'f lXM!, - (I - 2 Y)ci - 1] 


(34a) 


= - (1 - 2y)c| - Ij ■ 4>, - z®, <«*> 


(1 + m) 

'^ [Xlzllz-ll - 2y)c) - 1] 

- i 1 - 2 r) c ) - l] - 0 2 - z<t >2 


(34c) 

(34d) 


At this point a singular solution given by eqs (32) and (34) has been obtained which satisfies the main 
crack boundary conditions as well as the conditions at infinity. It is now necessary to model the branch. 


3.4. Branch Crack Problem 

The final potentials that give the solution to the interaction between a discrete dislocation and an 
interface crack subjected to uniform far-field stresses are: 
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(35a) 


- - D .03 - C 

<!>! = Oj + G>1 + Oi 

'F! = (35b) 

where (he Dislocation (D) potentials are given by equation (10), the infinity (°o) potentials are given by 
equation (20), and the Crack (C) potentials are given by equation (22). For the no interface case (a = P = 0) 
the present solution reduces analytically to Lo‘s (Lo 1978) solution. 

Replacing the discrete dislocation with a distribution of continuous dislocations and requiring this 
stress to be zero on the branch crack leads to an integral in terms of the unknown dislocation density. After 
non-dimensionalizing, separating the kernel into singular and regular parts, and some algebra, the following 
Cauchy-type integral equation is arrived at: 

f r jt s dt + l 4f ^ l(s ’ t)D(t) + ^ 2(8, t)0(t) ^ dt = _f(s) (36) 

in which D(t) = \i\ e*® 5/3t{ [u r ] + i [vq] }/irc(iq + 1) along the line zq = 1 + te 10 , [u r ] and [vq] are the 

jumps in the tangential and normal displacements across the dislocation line, Kj and K2 are the kernels 

resulting from the potentials derived earlier, and f(s) is the known traction along the branch line z=l + se* e 
due to the main crack loaded at infinity . 

This equation is solved numerically using piecewise continuous polynomials in the manner of 
Gerasoulis (1982). Once the dislocation densities have been determined numerically, the Stress Intensity 
Factors are obtained directly from the dislocation density at the tip of the branch crack in the normal 
fashion, e.g. Bryant, Miller, & Keer (1984). In applying the numerical scheme, the integration point at the 
base of the branch crack is ignored. Ignoring the integration point at the "knee" of the branched crack is an 
approximate method of incorporating the fact that the singularity at this end of the crack Is less than 1/2. 

He and Hutchinson (1988) investigated the effect of neglecting this integration point as opposed to 
including the actual singularity explicitly and found the effect to be minimal. A finite element analysis was 
run as part of the present study to check that the crack knee is indeed open under tensile loadings, and this 
was found to be the case for the crack geometries considered here. Thus, despite (lie approximate nature of 
the handling of the junction of the main crack with the branch, the numerical solution provides accurate 
results for the branching problem. 


4. Numerical Results and Conclusions 

The integral equation (36) was solved numerically for various crack branch geometries and material 
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property combinations for a pure tensile load at infinity, o yy . Following He and Hutchinson (1988), the 

parameter of interest is taken as g/go'- the ratio of the energy release rate at the tip of the branch to the 
energy release rate at an equivalent unbranched interface crack. For this analysis, Q corresponds to the strain 
energy release rate of a branch crack tip in material 1 and g$ corresponds to the strain energy release rate at 
the tip of an interface crack of total length 2c 1 = 2c + /. The expressions for Q , and g/go are: 


Q = 


1 - vi 

2 Hi 


(K? + K&) 


(37a) 
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(37b) 


Go 


L \ V JL + 1 - v 2 I (°»y + <&)(! + 4 £ 2 )tc c ' 
Hi J 4cosh 2 (7i£) 


£» 


(1 + a) 


cosh 2 7i e 


(l + 4e : 


: )k 2 


V 


Ki 

in c 


K u 

inc 


(37c) 


Figures 2 and 3 show the dependence of the strain energy release rate, gig o, on the branch angle, 9, for 
the cases of growth into the relatively stiff material and into the relatively soft material, respectively, for 
a = 0.5, P = 0, and 1/c = 0.001. As the length of the branch becomes very small relative to the main crack, 
it is expected that the present results should in some sense approach those of the semi-infinite crack. The 
gig o results in Figures 2 and 3 indeed agree well with those of He and Hutchinson (1988), although in 

Figure 2 g!g$ goes to 1 near 0 = 0 more quickly. 

For p ^ 0, there are some interesting problems regarding the loading parameters. As shown by 
Rice (1988), the far field loading and crack length are coupled in the expression for the complex intensity 

factor K = K\ + iK 2 = (kj + ik 2 ) V n cosh(rte), where k| + ik 2 is the complex intensity factor of Sih 
and Rice ( 1964). Since He and Hutchinson (1988) used y = tan* 1 KjfKy as their loading factor, it is not 
possible to match their loading parameter uniquely for p * 0. Specifically, the loading parameter used here, 

. oo <so 

y = tan 1 x y */a yy , is related to He and Hutchinson's parameter, y, by: 


= Y + tan 1 2 e - £ log 2 c (38) 

In the present analysis, the salient length parameter is the ratio of the branch crack length to the main crack 
length, //c, and the absolute value of c is arbitrary. Since the value of c is arbitrary, \|f and y cannot be 
related in a definitive manner. Equation (38) shows that for P = 0, £ will be 0 and y is equivalent to y, and 
thus the results in Figures 2 and 3 can be directly compared to He and Hutchinson's. Conversely for p ^ 0, 
y it y ; the two loading parameters are different. As noted by Rice ( 1988), £ is usually very small. 
However, because of the £ log 2c term in equation (38), this does not ensure that \y = y. If £ log 2c « 1, 
then y s y. 


Comparing Figure 2 and Figure 3 shows that for the case of the crack growing into the softer material, 
the maximum £/£b, g/gomw occurs at an angle somewhat off the interface and in the softer material 

(0 > 0). Although it is not shown by these figures, the angle at which K\ is maximized in the softer 
material corresponds to a nearly zero value for Kn. These results are typical for different //c s and Ts, and 
are consistent with the observations of He and Hutchinson concerning the tendency of branching to occur in 
the softer material provided the material toughnesses are comparable. Subsequent results are presented for 0 
in the soft material only. 

Figure 4 shows gig b vs. 0 for various T and //c = 0.001. Notice that as T increases, so does g/g 0 max . 
The angle that maximizes glgo, 0 ma x» also increases as T increases. Again, both these trends agree with He 
and Hutchinson's (1988) results for y # 0, although for the reasons given above, direct comparison is no 
longer possible. It should be noted that the in-plane stresses, a**, are increased also by increasing material 
mismatch, although these stresses appear to have only secondary influence on the branch behavior. 
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Figures 5 and 6 show the influence of finite main crack size by considering (j/cjo vs. 0 for various // c 
values. Observe that in both Figures 5 and 6 as // c increases, g/gomax decreases towards a value of unity, 
while the branch angle corresponding to maximum energy release rate decreases. (Note, however, that in 
calculating g/g q, go is adjusted to be the strain energy release rate for a crack of length 2c + / . Thus 
<5ma* itself is not necessarily decreasing with increasing branch length, since g$ is increasing with /. In 
fact, g maLX reaches a minimum and then starts to increase again as // c is increased). Using g/gomax as a 

crack growth criterion, the decreasing 0 m ax with increasing branch length can be taken as an indication that 
the branch would tend to turn back parallel to the interface as it grows, although this is a crude way to 
predict crack trajectories (and can be shown to be incorrect for some cases: see Rubinstein, 1989), . Given 
the further result that the maximum energy release rate at the branch tip approaches that of the interface 
crack itself as the branch grows, then together this would imply that in a limited sense branching Ls 
irrelevant, since the driving force would be similar whether branching occurs or not. The critical issue, 
however, remains the relative toughnesses of the materials involved 

The primary conclusions to be drawn from the present analysis are that branching from an interface is a 
viable mode of crack growth, and that such branching is likely to occur initially at an angle somewhere 
between 10° and 47° in the softer of the two materials, provided the material toughnesses are similar. As 
the branch extends its growth characteristics will be altered, resulting in a tendency to return to a path 
parallel to the interface, with a driving force similar to that of an unbranched crack. 
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Figure 4 - Strain energy release rate, g/Gg versus extension angle, 0 , for 
//c=0.001, Vj=v 2 =0.3, r=l, 3, 10, 100 (a = 0, 0.5, 0.81818, 
0.98020; p = 0, 0.14286, 0.23377, 0.28006). 
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Figure 5 - Strain energy release rate, versus extension angle, 0 , for 
v, =v 2 = 0.3, T= 3, //c=0.001, 0.01, 0.1, 0.5, 1.0 (a = 0.5, 
p = 0.14286). 
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Figure 6 - Strain energy release rate, g/<^ versus extension angle, 0 , for 

v, =v 2 = 0.3, T= 10, //c=0.001, 0.01, 0.1, 0.5, 1.0 (a = 0.81818, 
|3 = 0.23377). 


21 



NASA 

National Aeronautics and 
Space Administration 

1. Report No. 

NASA CR- 1 82273 

4. Title and Subtitle 


Analysis of Interface Crack Branching 


Report Documentation Page 


2. Government Accession No. 


3. Recipient’s Catalog No. 


5. Report Date 

March 1989 

6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 

505-63-31 

9. Performing Organization Name and Address 

11. Contract or Grant No. 

Case Western Reserve University 

Department of Civil Engineering NAG3-856 

Cleveland, Ohio 44106 13 . jype 0 f Report and Period Covered 

— Contractor Report 

12. Sponsoring Agency Name and Address Final 

National Aeronautics and Space Administration -J4 Sponsoring Agency Code 

Lewis Research Center 

Cleveland, Ohio 44135-3191 

15. Supplementary Notes 

Project Manager, Bernard Gross, Structures Division, NASA Lewis Research Center. R. Ballarini, Case Western 
Reserve University; D.J. Mukai and G.R. Miller, Department of Civil Engineering, University of Washington, 
FX10, Seattle, Washington 98195. 

16. Abstract 

A solution is presented for the problem of a finite-length crack branching off the interface between two bonded 
dissimilar isotropic materials. Results are presented in terms of the ratio of the energy release rate of a branched 
interface crack to the energy release rate of a straight interface crack with the same total length. It is found that 
this ratio reaches a maximum when the interface crack branches into the softer material. Longer branches tend to 
have smaller maximum energy release rate ratio angles indicating that all else being equal, a branch crack will 
tend to turn back parallel to the interface as it grows. 


7. Author(s) 

R. Ballarini, D.J. Mukai, and G.R. Miller 


17. Key Words (Suggested by Author(s)) 

Crack branching 
Fracture mechanics 
Interface crack 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 39 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No of pages 

23 


22. Price* 


NASA FOHM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 




